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一 种 结合 像素 空间 信息 和 多 维 直方 图 的 
彩色 图 像 快速 分 割 算法 * 
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(西安 交通 大 学 理学 院 ， 西 安 710049) 


H5 要 : 本 文 提 出 了 一 种 新 的 结合 空间 和 彩色 信息 的 彩色 图 像 快 速 分 割 方法 ， 首 先 对 三 个 单 色 的 直方 图 进 
行 多 阔 值 分 割 得 到 一 组 谷 值 和 峰值 ， 再 用 谷 值 向 量 将 真 彩色 三 维 直 方 图 分 割 成 不 同 区 域 ， 然 后 以 
峰值 组 成 种 子 向 量 实现 三 维 直 方 图 的 多 阔 值 分 割 ， 最 后 基于 空间 相 邻 关系 的 最 优 准则 对 分 割 的 区 
域 不 连续 性 和 “ 洞 ” 重 新 分 割 ， 提 高 了 过 渡 区 域 像素 的 分 割 精度 . 

关键 词 : 彩色 图 像 分 割 ， 真 彩色 三 维 直 方 图 ， 均值 漂移 算法 ， 多 阔 值 分 割 
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彩色 图 像 分 割 就 是 利用 图 像 的 彩色 和 空间 信息 ， 将 图 像 分 割 成 一 些 感 兴趣 的 区 域 . 灰 度 
图 像 分 割 方法 比较 成 熟 " 引 ， 但 将 其 直接 用 于 三 维特 征 空间 且 内 部 聚 类 不 紧凑 的 彩色 图 像 
分 割 ， 结 果 很 不 理想 . Ong 和 Hew 加 提出 了 和 迭代 阔 值 分 割 算 法 ， 他 们 在 HSI 彩 色 空 间 中 ， 
首先 根据 灰 度 直方 图 对 3 个 分 量 进行 阐 值 分 割 ， 求 出 各 分 割 区 域 中 象 素 的 RGB 分 量 的 均 
值 ， 并 进一步 求 出 相 令 区域 间 的 颜色 距离 ， 然 后 根据 给 定 的 门槛 值 进行 区 域 融 合 ， 以 消除 
碎片 Sun 和 Horng 中 提出 的 自动 选择 阐 值 的 彩色 图 像 分 割 算 法 只 适合 对 同一 类 的 图 像 进 
行 处 理 ， 且 易 受 环境 因素 的 影响 ，Kungollugb 等 将 真 彩 色 三 维 直 方 图 投影 成 三 个 二 维 直方 
图 (RG,GB,BR)， 利 用 三 个 二 维 直方 图 分 别 对 图 像 进行 多 阔 值 分 制 ， 再 融合 分 割 图 像 得 到 最 
终 分 割 结果 . 潘 晨 和 顾 锋 回 利用 尺度 空间 聚 类 理论 和 区 域 生长 法 对 二 维 多 阔 值 方法 进行 修改 ， 
降低 了 计算 复杂 度 ， 但 在 实现 过 程 中 仍然 需要 人 工 于 预 ， 调 整 分 割 维 数 、 分 辨 率 和 腐蚀 膨胀 次 
数 ， 叶 齐 祥 等 6 用 颜色 粗糙 度 对 图 像 进行 颜色 量化 后 使 用 增 量 式 的 区 域 生长 算法 形成 图 像 的 初 
始 分 割 区 域 ， 再 对 初始 分 割 区 域 进行 分 级 合并 ， 最 后 利用 形态 学 的 有 关 算 法 对 分 割 区 域 的 边缘 
进行 平滑 . 但 这 两 种 方法 分 割 图 像 的 后 处 理 比 较 复 杂 . 本文 利用 多 阐 值 算法 对 三 个 单 色 直方 图 
得 到 各 自 的 谷 值 和 峰值 ， 由 谷 值 组 合 分 割 真 彩色 的 三 维 直方 图 得 到 一 组 分 割 区 域 ， 以 峰值 组 合 
作为 相应 区 域 的 种 子 点 ， 在 真 彩色 的 三 维 直方 图 上 应 用 基于 均值 漂移 (Mean Shift, MS) 算法 得 
到 三 维 直方 图 多 立 值 分 割 的 类 中 心 颜色 向 量 ， 并 对 图 像 像 素 进行 类 标记 得 到 类 图 ， 然 后 对 类 图 
依据 最 优 分 割 准则 ， 结 合 空间 信息 得 到 改进 后 的 分 割 图 像 . 本 文 算法 应 用 血液 细胞 图 像 的 分 割 
试验 ， 结 果 表明 算法 的 有 效 性 . 
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2 ” 真 彩色 图 像 的 三 维 直 方 图 (3D Histogram) 的 分 割 


2.4 ”3D 直方 图 的 特性 

真 彩色 图 像 的 3D 直方 图 是 一 种 可 视 化 数据 图 ， 从 中 可 以 观察 到 像素 在 3D 颜色 空间 中 的 聚 
类 特性 . 利用 3D 直方 图 的 峰值 可 以 得 到 像素 分 类 的 各 类 中 心 的 颜色 向 量 ， 运 用 峰 态 间 的 谷底 
轮廓 曲线 形成 类 间 划 分 . 但 真 彩 色 的 3D 直方 图 存在 像素 数 分 布 密度 底 ， 聚 类 内 部 紧凑 性 差 ， 
分 布 空间 广 等 特点 ， 这 使 得 3D 直方 图 的 分 割 尤 为 困难 . 如 像素 在 空间 直 柱 上 的 分 布 可 能 会 集 
中 在 某 一 区 域 ， 从 而 存在 大 量 的 稀 疏 甚至 是 没有 像素 点 的 区 域 ， 例 如 256 x 256 大 小 的 图 像 共 
有 像素 65536 个 ， 而 3D 直方 图 的 颜色 柱 数 为 256 x 256 x 256 = 16777216， 绝 大 多 数 区 域内 的 
像素 分 布 稀少 . 用 一 般 的 优化 方法 寻找 直方 图 的 局 部 谷底 和 峰值 的 搜索 空间 太 大 ， 计 算 复杂 
高 ， 即 使 用 在 高 维 空间 比较 适用 的 均值 漂移 (MS) 算法 也 由 于 初始 种 子 点 的 选择 和 可 能 落 入 稀 
玻 区 域 而 使 算法 收敛 太 慢 ， 甚 至 不 收敛 . 对 3D 直方 图 的 多 阔 值 分 制 有 两 种 常见 的 方法 : 一 是 
用 三 个 单 色 直方 图 的 谷底 (有 图 maR 个 ，G 图 mc 个 ，B 图 ms 个 ) 形 成 一 组 平行 于 坐标 平面 的 
平面 艇 分 割 3D 直方 图 为 (ma 十 1) x (ma 9-1) x (mg 十 1) 块 区 域 ( 类 )， 并 以 块 内 像素 的 颜色 向 
量 均 值 为 类 中 心 的 颜色 向 量 ， 该 类 方法 使 图 像 过 度 分 割 ， 且 形成 了 的 类 中 心 颜色 可 能 在 图 像 上 
不 存在 .二 是 将 3D 直方 图 投影 到 三 个 坐标 平面 形成 三 个 二 维 直方 图 (RG,RB,GB)， 分 别 对 三 
个 二 维 直方 图 分 割 ， 然 后 融合 分 割 结果 得 到 3D 直方 图 的 分 割 ， 融 合 结果 取决 于 融合 准则 和 图 
像 数 据 ， 容 易 形 成 大 量 无 法 指派 的 像素 . 

2.0 MS 算法 的 缺陷 

MS 算法 是 一 种 利用 核 函 数 的 密度 估计 过 程 、 无 需 监督 和 参数 估计 ， 是 一 种 简单 有 效 的 聚 类 
方法 ， 有 自 适 应 的 步 长 ， 平 滑 的 搜索 轨迹 由， 最 速 上 升 的 性 质 ， 能 适合 估计 特征 空间 的 任意 形 
状 的 多 峰 聚 类 . 搜索 窗 宽度 参数 h 较 大 时 ，MS 征 空 间 中 密度 较 小 的 特征 具有 拟 制 作用 ， 因 此 
有 具有 较 少 特征 支持 的 类 可 能 被 忽略 而 收敛 于 其 附近 的 密度 较 大 的 峰 ， 对 于 特征 接近 的 聚 类 ， 也 
可 能 会 因 类 间距 (特征 中 心 向 量 的 距离 ) 小 于 搜索 窗口 宽度 而 被 MS 算法 归于 同一 类 . 当 户 较 小 
时 ，MS 对 特征 空间 的 局 部 聚 类 峰 比 较 敏 感 ， 容 易 受 噪声 的 干扰 . 均值 移动 向 量 的 计算 公式 确 
定 了 聚 类 的 时 效 性 . 假设 待 分 割 的 图 像素 数 为 N?， 每 个 像素 所 对 应 的 特征 问 量 的 维 数 为 4d， 则 
每 次 均值 移动 向 量 的 计算 需要 做 d x N? 比较 ， 若 每 个 点 的 平均 迭代 步 数 为 m， 则 图 像 分 制 所 
需 的 计算 量 为 m x d x N? x N2， 即 像素 数 的 平方 级 ， 对 于 稀疏 和 平缓 区 域 ，MS 算法 的 收敛 性 
很 差 ， 因 而 对 于 三 维 真 彩色 直方 图 难以 有 效 分 制 . MS 算法 虽然 有 这 些 缺 陷 ， 但 对 于 多 维 空间 
的 局 部 寻 优 具有 较 高 的 适用 性 和 重 棒 性 ， 为 此 本 文 对 MS 算法 进行 改进 ， 提 高 其 搜索 有 效 性 ， 
利用 其 寻找 细胞 彩色 图 像 3D 直方 图 的 分 割 阔 值 . 

2.3 ”基于 均值 漂移 的 3D 直方 图 多 模 态 峰 值 寻 优 

为 了 克服 MS 算法 在 平坦 区 域 收敛 绥 慢 和 模 态 峰值 收敛 点 依赖 种 子 点 的 缺陷 ， 我 们 把 基 
于 单 色 直 方 图 的 谷底 超 平 面 分 制 、 峰 值 组 合 向 量 和 MS 算法 结合 ， 提 出 一 种 基于 超 平面 分 割 
的 MS 算法 寻找 3D 直方 图 的 局 部 峰值 . 具体 做 法 是 : 

1) 先 对 图 像 的 三 个 单 色 直方 图 曲线 优化 ， 得 到 相应 的 谷 值 basini; 和 峰值 peakij， 其 
Pj = l,m k=l oeno m, Ri 为 第 i 个 彩色 直方 图 的 谷 值 个 数 和 峰值 个 数 ，i = 
1,2,3 表 示 红 、 绿 、 蓝 三 色 ; 

2) 由 x = basin, y = basing, z = basing (j = 1l,:,mi, k = l,,mg,l 一 
1,… ,m3) 2298 3D 直方 图 得 到 (ma +1) x (ma +1) x (ms 十 1) 个 空间 区 域 ; 

3) Hipeak,;, peakak, peaksl (j — 1 mi, k= 1, ng, 7 一 1 ,ns) 在 对 应 的 区 域 形 
成 (ni x no x n3) 个 种 子 点 ; 

4) 用 MS 算法 对 3) 形 成 的 所 有 种 子 点 进行 均值 漂移 ， 对 于 第 i 个 种 子 点 如 果 移 出 相应 区 域 
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沙 入 第 ;个 区 域 ， 则 合并 i,j 两 个 区 域 ， 并 从 种 子 点 集中 删 去 第 ?个 种 子 点 . 重复 步 又 4 直到 所 
有 种 子 点 收敛 ， 

通过 以 上 步骤 得 到 的 种 子 点 的 颜色 向 量 即 为 聚 类 的 类 中 心 ， 种 子 点 所 对 应 的 区 域 所 有 像素 
即 为 一 类 .为 克服 过 分 割 现象 ， 利 用 该 直方 图 空间 中 种 子 点 坐标 处 的 像素 频 度 值 ， 以 频 度 相同 
或 相近 为 连通 性 条 件 ， 在 原 3D 直方 图 空间 进行 三 维 区 域 生 长 ， 可 以 找 回 被 不 当 闭 值 切割 而 丢 
失 的 部 分 像素 ， 以 改善 分 割 结果 .由 于 直方 图 的 多 阔 值 分 割 没 有 考虑 图 像 的 空间 结构 ， 会 造成 
图 像 的 过 度 分 割 ， 因 而 需要 一 定 的 后 处 理 . 


3 ”基于 图 像 空 间 相 邻 关 系 和 最 优 准 则 的 再 分 割 


为 克服 3D 直方 图 阀 值 方法 只 考虑 彩色 信息 ， 而 不 考虑 各 像素 空间 信息 的 固有 缺点 提出 一 种 
基于 空间 相 邻 关系 的 最 优 分 割 准 则 ， 并 基于 此 准则 完成 相 邻 区 域 合 并 和 不 同类 别 区 域 过 渡 像 素 
的 划分 ， 消 除 阔 值 分 割 所 造成 的 目标 中 小 空洞 和 产生 的 其 他 小 碎片 . BU ELA EE Js ABEL TR DS 
个 过 程 : 空间 相 邻 区 域 的 合并 和 过 渡 区 像素 的 划分 . 

3.1 最 优 分 割 准 则 

图 像 I 的 像素 个 数 为 N, vz 为 像素 z = (x,y) € 工 的 颜色 向 量 (r,9， d), m 为 图 像 颜 色 向 
量 的 均值 . 图像 经 上 述 方法 分 割 后 成 为 几 个 大 类 ， 并 依据 图 像 像 素 的 空间 位 置 (二 维 坐 
bk) 分 成 M 个 区 域 ， 并 对 每 个 区 域 的 像素 赋予 不 同 的 类 标号 . Zi mi Ni 分 别 表示 第 i(i = 
1 .… ,AM) 个 区 域 的 像素 集 、 颜 色 均值 向 量 、 像 素数 .定义 如 下 公式 


1 
m=7 vs (1) 
zcl 
-xX (2) 
Mi N; Uz; 
z€Zi 
1 2 
$, = xw 2 les — ml (3) 
zcl 
1 
Si = N. » vz — mill , (4) 
k z€Zi 
M. 
Se=} Si (5) 
i—l 


其 中 5, 表示 图 像 的 平均 距离 ，S; 表示 第 ;个 区 域内 的 平均 距离 ，5- 表示 图 像 类 内 平均 距离 . 
定义 图 像 的 分 割 指标 为 
J = (S, m Se)/Se, (6) 
定义 图 像 某 个 局 部 区 域 的 分 割 指标 为 
Jy = (Sgk as Sa)/ Sek- (7) 


图 像 分 割 指标 了 反映 了 图 像 的 分 割 效 果 ， 图 像 分 割 效 果 好 ， 则 图 像 分 割 的 每 个 区 域内 ， 颜 
色 分 布 比较 均匀 ， 区 域内 平均 距离 就 校 则 指标 7 越 大 ， 否 则 指标 J 越 小 .为 了 反映 图 像 整 体 的 
分 割 效果 给 出 分 割 评价 指标 的 定义 
M 


1 
Jopt — N 2: 9 (8) 
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通过 最 小 化 公式 (8)， 可 以 实现 过 渡 区 域 像素 的 划分 . 但 是 这 一 优化 过 程 比较 复杂 ， 为 此 我 
们 采用 如 下 两 个 阶段 完成 . 

3.2 ”基于 分 割 最 优 准 则 的 空间 相 邻 区 域 合 并 

图 像 直方 图 阔 值 分 割 后 ， 并 依据 图 像 像 素 的 空间 位 置 (二 维 坐 标 ) 分 成 M 个 区 域 ， 并 对 每 个 
区 域 的 像素 赋予 不 同 的 类 标号 ， 把 区 域 看 成 一 个 点 ， 建 立 图 像 的 区 域 分 割 图 即 “ 类 -图 ”. 根 
据 区 域 的 空间 邻接 和 分 割 最 优 准 则 确定 两 个 相 邻 区 域 的 合并 ， 有 具体 步骤 如 下 : 

1) 建立 图 像 “ 类 -图 ”和 各 区 域 的 邻接 关系 表 ， 令 GZ 表 示 所 有 类 的 集合 ， 按 公式 (8) 计 
算 Jopt; 

2) EAZ, WAR. BWARRE ZP PREX, Z= 2 一 { 讨 ， 对 于 ?从 2 中 取出 与 ?有 邻 
接 关 系 的 类 组 成 类 i 的 邻接 集 Adj(i); 

3) AGOA, 5&2) 否则 ， 从 Adj( 让 中 取出 元 素 k，Adj(i) = Adj(i) — (kp. it 
& Z' = ZiU Zn， 对 集合 2Z' 按 公式 (8) 计算 新 的 Jopti; 

4) E Jop > Jopt1， 则 合并 类 i 和 类 k， 令 Jopts， 修 改 个 区 域 的 邻接 关系 ， 转 2); 和 否则 ， 转 
步骤 3)， 

经 过 上 述 步 又 完成 了 空间 相 邻 区 域 的 合并 ， 消 除了 分 割 造成 的 润 和 小 碎片 . 

3.3 ”类 间 过 度 区 域 像素 的 划分 

对 每 个 区 域 的 边界 内 2 像素 标记 为 过 度 像素 . 对 每 个 度 像 素 i 以 其 为 中 心 建立 一 个 半径 
为 7 的 邻 阀 ， 在 邻 域 中 然后 对 每 一 像素 作 松 弛 操作 ， 考 察 其 相 邻 的 区 域 ， 对 每 个 相 邻 区 域 把 像 
素 i 并 入 其 中 ， 计 算 该 邻 域 的 分 割 指标 Js， 比较 所 有 分 割 的 局 部 邻 域 指标 值 ， 把 最 小 指标 值 
对 应 的 像素 划分 为 像素 i 的 最 终 划 分 . 





4 ”融合 空间 信息 的 直方 图 多 阅 值 彩色 图 像 分 割 方法 及 步 又 


真 彩色 细胞 图 像 分 割 的 步骤 如 下 : 

1) 对 彩色 图 像 生成 真 彩色 三 维 直方 图 和 三 个 单 色 直 方 图 ， 并 对 直方 图 进行 光滑 ; 

2) 对 三 个 单 色 直方 图 寻求 谷底 和 峰值 

3) 以 单 色 直方 图 的 谷底 值 对 三 维 直 方 图 进行 分 割 ， 以 单 色 直方 图 的 峰值 组 合成 三 维 空间 
的 种 子 向 量 ， 并 与 前 述 的 分 割 区 域 相对 应 ; 

4) 按照 2.3 所 述 的 方法 基于 均值 漂移 进行 各 种 子 点 的 局 部 峰值 寻 优 ， 得 到 各 区 域 的 聚 类 中 
心 颜色 向 量 ; 

5) ” 按 区 域 类 中 心 对 各 区 域 标号 生成 图 像 “ 类 -图 ”; 

6) 计算 图 像 的 分 割 指标 值 ， 并 3.2 方 法 进行 区 域 合并 ; 

7) ”依据 3.3 对 过 度 像素 进行 划分 . 


5 ”细胞 彩色 图 像 的 分 割 步骤 及 实验 结果 


本 文 图 像 从 瑞 氏 染色 的 人 体 血液 涂 片 中 获得 . 图 1(a) 为 细胞 的 真 彩色 图 像 ， 图 1(b), 图 1(c)， 
图 1(q) 分别 为 对 应 的 单 色 直方 图 多 立 值 分 割 图 像 .在 原始 单 色 直方 图 中 寻找 合理 的 谷底 作为 分 
制 阐 值 ， 分 割 效 果 很 不 理想 ， 图 1(e) 是 用 本 文 提出 基于 均值 漂移 的 3D 直方 图 多 模 态 峰值 寻 优 
方法 分 割 成 的 七 类 图 像 ， 图 2(a) 是 在 分 割 基础 上 对 过 度 像 素 划 分 后 的 图 像 和 分 割 准则 指标 的 计 
算 结果 Jopt = 0.203; 图 2(b) 是 基于 最 优 分 割 准则 分 割 后 并 进行 区 域 合并 和 去 除 最 小 区 域 后 得 
到 的 分 割 图 像 ， 其 分 割 准则 指标 JS. = 0.158. 细胞 图 像 被 分 为 有 核 细胞 (细胞核 和 细胞 质 )， 
无 核 细 胞 和 背景 共 四 类 ， 且 各 细胞 图 像 连 成 一 个 整体 ， 基 本 满足 了 医学 上 对 细胞 识别 的 分 割 要 
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求 . 对 于 实验 图 像 ， 本 文 方法 分 割 成 4 个 有 核 细 胞 ，33 个 无 核 细 胞 ， 人 工 方法 分 割 成 4 个 有 核 
细胞 ，31 个 无 核 细胞 ， 对 血液 细胞 的 10 幅 图 像 用 本 文 方法 和 人 工 方法 进行 了 分 割 实验 ， 得 到 
了 的 统计 结果 表 1 (单位 为 细胞 的 个 数 ). 可 以 看 出 本 文 实验 结果 与 人 工 结果 基本 接近 ， 特 别 是 
有 核 细胞 的 分 割 结果 . 本文 实验 是 在 RGB 通道 上 进行 的 ， 而 彩色 转换 成 Lab 和 HIS 空间 ， 细 
胞 图 像 的 分 割 效 果 更 好 一 点 . 








(c) 绿色 分 割 图 像 


T 





(d) 蓝 色 分 割 图 像 (e) 基于 三 维 直方 图 的 多 模 态 分 害 
图 1: 细胞 彩色 图 像 的 分 割 图 像 
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(a) 分 割 图 像 过 度 像素 划分 后 (Jopt = 0.203) 


(b) 经 最 优 分 割 准则 的 区 域 合并 后 (Jopt = 0.158) 


图 2: 本 文 方法 分 割 的 细胞 图 像 


Aio 细胞 图 像 分 割 结果 比较 





















































图 像 编 号 1 2 3 4 5 6 | 7 8 | 9 | 10 
有 核 AL 2 E 5 3 6 8 3 7 5 6 
细胞 |xx|co2la4lsis|sls ESEL EDE 
无 核 | 人 工 | 21 |33| 41 | 30 | 65 | 90 | 42 | 45 | 47 | 72 
细胞 | 本 文 | 20 31 42 30 64 93 | 44 | 46 | 4v | 7O 

图 像 编号 11 12 13 14 15 16 17 18 19 | 20 
有 核 Ax | 7 5 4 3 3 TER 6 | N 2 
细胞 | AX 8 | 5 4 3 3 4 5 6 7 2 
无 核 人 工 78 | 66 | 43 | 32 | 36 | 41 | 32 | 56 | 63 | 23 
细胞 本 文 79 | 64 44 34 34 39 | 34 55 59 23 








本 文 提出 的 一 种 结合 像素 空间 信息 和 三 维 直 方 图 的 彩色 图 像 快 速 分 割 算法 . 首先 利用 多 阐 
值 算法 和 均值 漂移 算法 得 到 三 维 直 方 图 的 分 类 中 心 颜 色 向 量 ， 然 后 对 图 像 像素 进行 类 标记 得 
到 类 图 ， 最 后 对 类 图 依据 最 优 分 割 准则 并 结合 空间 信息 改进 分 割 图 像 . 算法 消除 了 分 割 中 的 
“ 洞 ” 和 细小 碎片 ， 提 高 了 过 渡 像 素 的 分 割 精度 ， 为 图 像 的 细胞 识别 莫 定 了 基础 . 
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A Fast Color Image Segmentation Algorithm by Combining Pixel 
Space Information and 3D Histogram 


ZHANG Gai-ying, ZHANG Jiang-she 


(School of Science, Xi'an Jiaotong University, Xi'an 710049) 


Abstract: By combining the pixel space information and color information, a fast color image seg- 
mentation algorithm is proposed in this paper to segment blood cell color images. The main ideas 
of this method are as follows: firstly, we apply the multi-threshold segmentation method to the three 
monochrome histograms to obtain a set of combinations of valley and peak values. Then, we utilize 
valley values to construct some hyperplanes in order to segment the true-color 3D histogram into differ- 
ent areas and employ the mean-shift algorithm with the seed vectors that are composed of peak values 
to implement the multi-threshold segmentation of 3D histograms. In order to solve the discontinuity 
problem and “hole” phenomenon of the areas obtained by single-color segmentation, a re-segmentation 
method using the optimal criterion based on adjacent space information is suggested. Some exper- 
iments are conducted and the results show that it can improve the segmentation accuracy to some 
extent. 

Keywords: color image segmentation; true-color 3D histogram; mean-shift algorithm; multi-threshold 


segmentation 
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